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ABSTRACT 

The  monthly  mean  TOPEX/POSEIDON  crossover  data  in  the  South  China  Sea  are  used  to  investigate  the 
spatial  and  temporal  variability  of  sea  surface  height  anomaly  (SSHA).  Multi-time  scale  variability  is  found  using  the 
empirical  orthogonal  function  (EOF)  analysis.  The  seasonal  variability  dominates  SSHA  with  two  spatial  patterns: 
basin-wide  gyre  and  north-south  double  gyres.  The  intraseasonal  variability  of  SSH  has  high  spatial  variability  at  the 
continental  shelf  such  as  west  of  Hainan  Island  and  Borneo,  and  east-west  double  gyre  structure  in  the  deep  basin.  The 
interannual  variability  of  SSHA  has  a  north-south  double  gyre  pattern.  The  scale  interaction  between  seasonal  and 
interannual  processes  may  be  taken  place  at  the  north-south  double  gyre  pattern. 


Keywords:  TOPEX/POSEIDON  altimetry.  South  China  Sea,  empirical  orthogonal  functions,  seasonal,  intraseasonal, 
interannual  variability 


1.  INTRODUCTION 

The  South  China  Sea  (SCS)  is  a  semi-enclosed  tropical  and  sub-tropical  sea  located  between  the  Asian  landmasses 
(including  Taiwan)  to  the  west  and  north,  the  Philippines  Islands  to  the  east  and  Malaysia  and  Borneo  to  the  south 
(Figure  1),  a  total  of  3.5xl06  km2  surface  area.  It  includes  the  shallow  Gulf  of  Thailand  and  connections  to  the  East 
China  Sea  (through  the  Taiwan  Strait),  the  Pacific  Ocean  (through  the  Luzon  Strait),  the  Sulu  Sea  (through  the  Balabac 
Strait),  the  Java  Sea  (through  the  Gasper  and  Karimata  Straits)  and  to  the  Indian  Ocean  (through  the  Strait  of  Malacca). 
All  of  these  straits  are  shallow  except  the  Luzon  Strait  whose  maximum  depth  is  1800  m.  Consequently  the  SCS  is 
considered  a  semi-enclosed  water  body.  The  complex  topography  includes  the  broad  shallows  of  the  Sunda  Shelf  in  the 
south  and  southwest;  the  continental  shelf  of  the  Asian  landmass  in  the  west  and  north,  extending  from  the  Gulf  of 
Tonkin  to  the  Taiwan  Strait;  a  deep,  elliptical  shaped  basin  in  the  center,  and  numerous  reef  islands  and  underwater 
plateaus  scattered  throughout.  The  shelf  that  extends  from  the  Gulf  of  Tonkin  to  the  Taiwan  Strait  is  consistently  about 
70  m  deep,  and  averages  150  km  in  width  from  the  shoreline;  the  central  deep  basin  is  1900  km  along  its  major  axis 
(northeast-southwest)  and  approximately  1100  km  along  its  minor  axis,  and  extends  to  over  4000  m  deep.  The  Sunda 
Shelf  is  the  submerged  connection  between  Southeast  Asia,  Malaysia,  Sumatra,  Java  and  Borneo  and  is  100  m  deep  in 
the  middle;  the  Gulf  of  Thailand  is  about  70  m  deep  in  its  central  part. 

The  observed  circulation  patterns  of  the  intermediate  to  upper  layers  of  the  SCS  are  also  primarily  forced  by  the  local 
monsoon  systems.  Based  on  limited  data,  Wyrtki  (1961)  constructed  the  surface  circulation  pattern.  From  April  to 
August,  the  weaker  southwesterly  summer  monsoon  winds  result  in  a  wind  stress  of  over  0.1  N/m2  (Figure  2a),  which 
drives  a  northward  coastal  jet  off  Vietnam  and  an  anticyclonic  circulation  in  the  central  SCS.  From  November  to  March, 
the  stronger  northeasterly  winter  monsoon  winds  corresponding  to  a  maximum  wind  stress  of  nearly  0.2  N/m2  (Figure 
2b)  causing  a  southward  coastal  jet  off  Vietnam  and  a  cyclonic  circulation  in  the  central  SCS.  Transitional  periods  are 
highly  variable  winds  and  surface  currents.  The  current  meter  observation  in  SCS  is  very  sparse.  It  is  very  hard  to  get 
spatial  and  temporal  variability  of  the  SCS  circulation  from  in-situ  data.  Sea  surface  elevation  anomaly  data  from  the 
TOPEX/POSEIDON  (T/P)  provide  useful  information  about  the  surface  circulation. 
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Figure  1.  Geography  and  isobaths  showing  the  bottom  topography  of  the  South  China  Sea. 
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Figure  2.  Climatological  wind  stress  for  (a)  June  and  (b)  December  (After  Hellerman  and  Rosenstein  1983). 


2.  TOPEX/POSEIDON 

T/P  is  a  joint  venture  by  NASA  and  the  French  space  Centre  National  d’Etues  Spatiales  (CNES)  that  is  designed  to 
produce  global  maps  of  ocean  surface  topography  (Stewart  et  al.,  1986).  The  data  used  in  this  study  are  from  October 
1992  to  September  1999  (257  cycles).  There  are  24  crossover  points  over  SCS  (Figure  3).  After  removing  the  tides,  the 
monthly  mean  values  (84  months  all  together)  are  calculated  for  each  crossover  point. 


3.  EMPIRICAL  ORTHOGONAL  LUNCTION  (EOL)  ANALYSIS 

3.1.  Principle  EOL  Modes 

The  monthly  sea  surface  height  anomaly  (SSHA)  data  from  T/P  form  an  NxP  matrix.  Here  P  =  84,  is  the  total  number  of 
time  points  used  for  computing  the  covariance  matrix.  N  =  24,  corresponds  to  the  number  of  crossover  points.  Empirical 
Orthogonal  Function  (EOF)  analysis  (Lorenz,  1956)  widely  used  in  oceanographic  and  meteorological  research  is  the 
same  as  Principal  Component  (PC)  Analysis  (Hotelling,  1933)  in  the  statistics  community.  PCs  are  the  amplitudes, 
which  are  functions  of  time,  of  their  corresponding  EOFs.  These  EOFs  can  be  found  by  calculating  the  unitary 
eigenvectors  of  the  covariance  matrix  associated  with  the  sample  data  field.  EOF  analysis  separates  the  data  sets  into 
eigenmodes.  Generally  speaking,  each  mode  has  an  associated  variance,  nondimensional  spatial  pattern,  and  dimensional 
time  series.  Therefore,  we  have  84  time  points  in  the  EOF  analysis.  From  this  data  matrix  a  24-square  spatial  covariance 
matrix  is  calculated. 

SSHA  can  be  quantitatively  investigated  with  the  method  of  EOF  analysis  and  a  number  of  its  generalized  forms.  In 
order  to  delineate  the  major  modes  of  variability,  we  perform  an  EOF  analysis  using  the  method  outlined  in  Section  3 
and  obtain  the  first  five  leading  EOFs,  which  are  able  to  account  for  almost  66.9%  of  the  total  variance  during  the 
seasons  (Table  1).  When  combined  EOF1,  EOF2  and  EOF3  explain  near  half  (49.5%)  of  the  total  variance.  Each  EOF 
mode  is  normalized  so  that  its  total  spatial  variance  is  equal  to  unity.  So,  those  patterns  of  the  first  five  EOFs  are  enough 
to  explain  the  spatial  anomalies  of  the  whole  SCS.  Hence,  it  suffices  to  focus  on  the  first  five  EOFs. 

3.2.  EOF1 

The  EOF1  mode  (Figure  4a)  shows  a  basin-scale  SSHA  structure  with  high  values  (larger  that  0.1)  in  the  SCS  deep  basin 
and  low  values  near  Chinese  and  Vietnamese  coast.  The  surface  circulation  corresponding  to  EOF1  is  a  basin-scale  gyre. 
To  identify  the  gyre  is  cyclonic  or  anticyclonic,  we  should  consider  the  product  of  the  EOF  values  and  the  corresponding 
time  series  of  the  amplitudes,  also  called  PC  components. 


The  first  principal  component,  PCi,  for  October  1992  -  September  1999  is  shown  in  Figure  4b.  The  time  series  of  PCi 
(Figure  4b)  shows  the  seasonal  variability.  The  Fourier  analysis  on  PCi  shows  the  existence  of  a  strong  peak  at  one 
year’s  period  (Figure  4c).  During  the  winter  monsoon  season,  PCi>0,  the  first  mode  shows  positive  SSF1A  in  the  SCS 
deep  basin,  which  implies  an  anticyclonic  gyre  in  SCS.  During  the  summer  monsoon  season,  PCi<0,  the  first  mode 
shows  negative  SSF1A  in  the  SCS  deep  basin,  which  implies  a  cyclonic  gyre. 

Table  1.  Variances  of  the  First  Five  Leading  EOFs 


EOF 

Variance  (%) 

Cumulative 
Variance  (%) 

1 

24.16 

24.16 

2 

14.59 

38.75 

3 

10.76 

49.51 

4 

6.77 

56.28 

5 

5.04 

61.32 

364  121 


Figure  4.  (a)  The  EOF1  mode  (unit:  cm),  (b)  time  series  of  PC  1?  and  (c)  Fourier  spectrum  of  PC]. 


3.3.  EOF2 


The  EOF2  mode  (Figure  5a)  shows  double-gyre  structure  with  a  high  center  in  the  southern  SCS  (south  of  10°N)  and  a 
low  center  in  the  northern  SCS  (north  of  10°N).  To  identify  each  gyre  is  cyclonic  or  anticyclonic,  we  should  consider  the 
product  of  the  EOF2  values  and  the  corresponding  time  series  of  PC2. 

The  second  principal  component,  PC2,  for  October  1992  -  September  1999  is  shown  in  Figure  5b.  The  Fourier  analysis 
on  PC2  shows  the  existence  of  a  strong  peak  (0.48)  at  one  year’s  period.  During  the  winter  monsoon  season,  PC2<0,  the 
second  mode  shows  low  SSFIA  center  in  the  southern  SCS  and  high  SSF1A  center  in  the  northern  SCS,  which  implies  a 
cyclonic  gyre  in  the  southern  SCS  and  an  anticyclonic  gyre  in  the  northern  SCS.  During  the  summer  monsoon  season, 
PC2>0,  the  second  mode  shows  high  SSFIA  center  in  the  southern  SCS  and  low  SSFIA  center  in  the  northern  SCS,  which 
implies  an  anticyclonic  gyre  in  the  southern  SCS  and  a  cyclonic  gyre  in  the  northern  SCS. 

Besides  the  seasonal  variability,  PC2  has  a  second  peak  (0.07)  at  849  days  (around  2  years  10  months).  Since  the  data 
length  is  7  years.  This  indicates  that  the  interannual  variability  contains  EOF2  pattern  (north-south  double  gyres). 


3.4.  EOF3 


The  EOF3  mode  (Figure  6a)  shows  weak  spatial  variability  (east-west  gyres)  in  the  SCS  deep  basin  and  high  spatial 
variability  at  the  continental  shelf  such  as  west  of  Flainan  Island  and  Borneo.  The  third  principal  component,  PC3,  for 
October  1992  -  September  1999  is  shown  in  Figure  6b.  The  Fourier  analysis  on  PC3  (Figure  6c)  shows  the  existence  of 
a  strong  peak  (0.08)  at  46days  (intraseasonal  time  scale).  Thus,  the  evident  intraseasonal  variability  of  SSHA  occurs  in 
the  continental  shelf  region. 


Figure  6.  (a)  The  EOF3  mode  (unit:  cm),  (b)  time  series  of  PC3,  and  (c)  Fourier  spectrum  of  PC 


3.5.  EOF4 


The  E0F4  mode  (Figure  7a)  shows  double-gyre  structure  with  a  low  center  in  the  southern  SCS  (south  of  10°N)  and  a 
high  center  in  the  northern  SCS  (north  of  10°N).  This  mode  is  similar  to  EOF2  mode.  The  time  series  of  PC4,  for  October 
1992  -  September  1999,  is  shown  in  Figure  7b.  The  Fourier  analysis  on  PC4  (Figure  7c)  shows  the  existence  of  a  strong 
peak  (0.15)  at  1019  days  (around  2  years  10  months).  Thus,  the  evident  interannual  variability  of  SSFIA  reveals  the 
north-south  double  gyre  structure.  The  second  peak  (0.1 1)  is  located  at  340  days  (seasonal  variability). 


3.6.  EOF5 


The  EOF5  mode  (Figure  8a)  shows  high  spatial  variability  at  the  continental  shelf  such  as  west  of  Flainan  Island  and 
Borneo,  and  east-west  double  gyre  structure  in  the  deep  basin.  The  fifth  principal  component,  PC5,  for  October  1992  - 
September  1999  is  shown  in  Figure  8b.  The  Fourier  analysis  on  PC,  (Figure  8c)  shows  the  existence  of  a  strong  peak 
(0.16)  at  85  days  (intraseasonal  time  scale). 


Figure  8.  (a)  The  EOF5  mode  (unit:  cm),  (b)  time  series  of  PC5,  and  (c)  Fourier  spectrum  of  PC?. 


4.  MULTI-SCALE  VARIABILITY 

T/P  data  show  the  existence  of  seasonal,  intraseasonal,  and  interannual  variabilities  of  SSHA. 

4.1.  Seasonal  Variability 

Seasonal  variability  dominates  SSHA  with  two  spatial  patterns:  basin-wide  gyre  (EOF1)  and  north-south  double  gyres 
(EOF2  and  EOF4).  The  seasonal  variability  of  basin-wide  gyre  is:  an  anticyclonic  gyre  during  the  winter  monsoon 
season,  and  a  cyclonic  gyre  during  the  summer  monsoon  season.  The  seasonal  variability  of  north-south  double  gyre  is: 


a  cyclonic  gyre  in  the  southern  SCS  and  an  anticyclonic  gyre  in  the  northern  SCS  during  winter  monsoon  season,  and 
an  anticyclonic  gyre  in  the  southern  SCS  and  a  cyclonic  gyre  in  the  northern  SCS  during  the  summer  monsoon  season. 

4.2.  Interannual  Variability 

Interannual  variability  of  SSHA  has  a  north-south  double  gyre  pattern  (EOF2  and  EOF4).  Since  this  pattern  experiences 
both  seasonal  and  interannual  variabilities,  the  scale  interaction  between  seasonal  and  interannual  processes  may  be 
taken  place  at  this  spatial  pattern,  which  is  the  out-of-phase  of  SSF1A  in  northern  and  southern  SCS  (increase  of  SSF1  in 
northern  SCS  and  decrease  of  SSFI  in  the  southern  SCS,  and  vice  versa). 

4.3.  Intraseasonal  Variability 

Intraseasonal  variability  of  SSFI  has  high  spatial  variability  at  the  continental  shelf  such  as  west  of  Flainan 
Island  and  Borneo,  and  east-west  double  gyre  structure  in  the  deep  basin  (EOF3  and  EOF5).  Shelf  processes  are  major 
sources  for  the  intraseasonal  variability. 


5.  CONCLUSIONS 

The  EOF  analysis  was  applied  to  the  monthly  mean  T/P  crossover  data  in  the  South  China  Sea  to  investigate  the  spatial 
and  temporal  variabilities  of  SSHA.  Our  study  shows  the  existence  of  multi-time  scale  variability:  seasonal,  interannual, 
and  intraseasonal.  Following  results  were  obtained. 

(1)  Seasonal  variability  dominates  SSHA  with  two  spatial  patterns:  basin-wide  gyre  and  north-south  double  gyres. 

(2)  Interannual  variability  of  SSHA  has  a  north-south  double  gyre  pattern. 

(3)  Intraseasonal  variability  of  SSH  has  high  spatial  variability  at  the  continental  shelf  such  as  west  of  Hainan  Island  and 
Borneo,  and  east-west  double  gyre  structure  in  the  deep  basin. 

(4)  Scale  interaction  between  seasonal  and  interannual  processes  may  be  taken  place  at  the  north-south  double  gyre 
pattern. 
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